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Abstract 

Consider a semi-infinite skew-symmetric moment matrix, evolv- 
ing according to the vector fields dm/dt^ = A k m + mA Tfc , where 
A is the shift matrix. Then The skew-Borel decomposition := 
Q~ 1 JQ T ~ 1 leads to the so-called Pfaff Lattice, which is integrable, by 
virtue of the AKS theorem, for a splitting involving the affine sym- 
plectic algebra. The tau-functions for the system are shown to be 
pfaffians and the wave vectors skew-orthogonal polynomials; we give 
their explicit form in terms of moments. This system plays an im- 
portant role in symmetric and symplectic matrix models and in the 
theory of random matrices (beta=l or 4). 
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Lie algebra splitting. Throughout this paper the Lie algebra T> = 
gloo of semi-infinite matrices is viewed as composed of 2 x 2 blocks. It 
admits the natural decomposition into subalgebras: 



(0.1) 



where T>q has 2x2 blocks along the diagonal with zeroes everywhere 
else and where T> + (resp. £>_) is the subalgebra of upper-triangular 
(resp. lower-triangular) matrices with 2x2 zero matrices along T>q 
and zero below (resp. above). As pointed out in (0.1), T>q can further 
be decomposed into two Lie subalgebras: 

Vq = {all 2 x 2 blocks € T>q are proportional to Id} 

£>+ = {all 2x2 blocks € T> have trace }. (0.2) 

Consider the semi-infinite skew-symmetric matrix J, zero everywhere, 
except for the following 2x2 blocks, along the "diagonal", 



/ 1 
-1 



J 



\ 



1 





1 





V 



€ Pq", with J 2 



-I, (0.3) 



and the associated Lie algebra order 2 involution 

J : V — >V:a* — ► J (a) := Ja T J. (0.4) 
The splitting into two Lie subalgebras[] 

V = k + n, with k = Z?_ +T>q and n = {a + J a, a £ T>} = sp(oo), 

(0.5) 



1 Note n is the fixed point set of J . 
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with corresponding Lie groups0<?k and Q n = Sp(oo), will play a crucial 
role here. Let ir^ and ir n be the projections onto k and n. Notice 
that n =sp(oo) and Q n = Sp(oo) stand for the infinite rank affine 



symplectic algebra and group; e.g. see [10]; the exposition in this 
paper is totally self-contained and does not depend on any knowledge 
of the affine symplectic algebra. 

Pfaff Lattice. As will be shown in this paper, the Lax pair, which 
we call the Pfaff lattice^], 

dL tr_Z7 +1 

= [-7r k V7^, L] = [n n V7U, L) , with Hi = — — , (0.6) 
oti i + 1 

(Pfaff Lattice) 

on matrices L = QAQ^ 1 , with Q £ Gk and A the customary shift 
operator, is completely integrable, as a result of the AKS-theorem. 
The Lax pair, written in compact form in (0.6), is given explicitly in 
(1.4). In pj], this lattice was investigated from the point of view of the 
2d-Toda lattice with special initial conditions. In Q, we developped 
the wave- and tau-function theory for this lattice and we exhibited a 
map from the Toda to the Pfaff lattice. 

Linearizing vector fields and Pfaffian f-functions. The key to 

this system is its linearization (moment map): namely, if L = QKQ~ l 
flows according to (0.6), then the following skew-symmetric matrix, 
constructed from Q, 

moo := Q- 1 JQ T - 1 

flows as 

^£ = A fc m 00 + m 00 A Tfc ,A: = 1,2,3,... . (0.7) 

Having constructed the skew-symmetric matrix m^it) = (fJ-ij(t)) Q<i j <OQ i 
we define the skew-symmetric submatrices 

m n (t) = (fHj(t)) Q<i j<n _x , with f 2n (t) = det 1/2 m 2 „(t), t G C°°. 

(0.8) 



2 {?k is the group of invertible elements in k, i.e., lower-triangular matrices, with non- 
zero 2x2 blocks proportional to Id along the diagonal. 

3 The Hamiltonians 7ii are viewed as formal sums; the convergence of this formal sum 
would require some sufficiently fast decay of the entries of L. Since V7i,; = L l , one does 
not need to be concerned about this point. 
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Note that, since m n is skew-symmetric, det m2 n +i = 0, whereas det rri2n 
is a perfect square, with det 1//2 m,2 n being the Pfaffian ; for a definition, 
specifying the sign, see section 3. The technology of letting rrioo flow 
in time and Borel decomposing has been used quite extensively in 
other situations, like the standard Toda lattice or the 2-Toda lattice 

1- 

The so-called Pfaffian f -function is not a KP r-function, but enjoys 
different bilinear identities and Hirota-type bilinear equations Q, 

{T2n(t ~ [u]),T 2n (t ~[v])} 

+ - V- 1 )(f 2n ( y t - [u])f 2n (t - H) - f 2n (t)f 2n (t - [u] - [«])) 

u 2 

= uv(u - v)f 2n -2(t - [u] - [v])f 2n+2 (t), where [u] = (u, ...). 
In particular, 

/ - 1 d 2 \ 

Pk+4{9) - TTTTT - ^ T 2n O T 2n = Pk{9) T 2n +2 ° T 2n -2, (0.9) 

^ I OtiOt k+3 J 

for k,n = 0,1,2,... . Both these equations are reminiscent of the dif- 
ferential Fay identities and the KP equations for usual r-functions 
O], but with an additional term on the right hand side, involving 
f 2n +2 and f 2n -2- 

For k = 0, this KP-like equation has already appeared in the con- 
text of the charged BKP hierarchy, studied by V. Kac and van de Leur 
plf l; the precise relationship between the charged BKP hierarchy of 
Kac and van de Leur and the Pfaff Lattice, introduced here, deserves 
further investigation. 

Skew-Borel decomposition and skew-orthogonal polynomials. 

Starting with equations (0.7), we show they have a unique semi-infinite 
skew-symmetric matrix solution moo(t), provided the initial condition 
satisfies detm^n 7^ 0, for n = 0,1,2,... It is given in terms of its 
(unique) "skew Borel decomposition"^] 

moo = Q~ 1 JQ T ~ 1 , with Q G Q k . (0.10) 

Here, we view moo as a matrix, providing the 11 inner-product" between 
monomials, 

{y l ,z j ) := (moo)ij =: /%•• (0.11) 



'unique, modulo a sequence of ±-signs. 
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Then the lower-triangular matrix Q is conveniently described in terms 
of a sequence of so-called " skew-orthonormal polynomials"^ 



<?0) = (qo(.z),qi{z),...) T = Qx{z) satisfying ((<?;, <?j))j 



M>0 



J. 



Skew-orthonormal polynomials were first introduced by Mehta [|T^]; 
see also Brezin and Neuberger ||, where they appeared in the context 
of unoriented random surfaces. To the best of our knowledge, neither 
the form, nor the connection with the Pfaff lattice was known. 

We show the polynomials q(t, z) in z, depending on t, are explicitly 
given by pfaffians of skew-symmetric matrices, which in the even case 
fen's) are formed by replacing in m-2 n +2 the 2n + 2th row and column 
by powers of z and in the odd case fen+i's), by replacing the 2n + 1th 
row and column by the same powers of z, keeping the skewness of the 
matrices: 



+2) 



-1/2 



Pf 



I o 

-Moi 

~M0,2n 

V -i 



Moi 




"Ml,2n 

—z 



Ml,2n Z 





~2n 



2n 

J 



q2n+l(t,z) 



{T2nT2n+2) 1/2 pf 



(T2nT2n+2) 1/2 ( Z + 7^- ) (f 2n T2n+2) i/2 <72n(*, Z). 



( ^01 

-M01 

-1 -z 

\— M0,2n+1 — Ml,2n+1 







\{z) = (l,z,z 2 ,...). 



1 M0,2n+1 \ 
^ Ml,2n+1 



-Z 2n+1 



2; 



2n+l 



/ 



(0.12) 
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Finally, the matrix 



L = QAQ , satisfying zq(t, z) = L(t)q(t, z), 



(0.13) 



forms a unique semi-infinite solution to the Lax pair (0.6). 

It is interesting to point out the striking similarity with the orthog- 
onal polynomials, the Toda lattice and the matrix models, described 
in H . Both, standard Toda[] and Pfaff lattices can be viewed as reduc- 
tions of the 2- Toda lattice, where the initial condition is a Hankel 
matrix for the standard Toda lattice and a skew-symmetric matrix for 
the Pfaff lattice. In the former case, the motion takes place within the 
big stratum and in the latter, within a deeper stratum. 



1 The Pfaff Lattice and a Lie algebra 
splitting 

Remember the decomposition (0.1) of the Lie algebra into subalgebras: 

v = gioo = z>- © £> e v+ = V- © Vq © v£ © v+. (1.1) 

Remember the element J € T>q , with J 2 = —I and the associated 



map: 



J : V — ► V : a i — ► J (a) := Ja T J. 



Theorem 1.1 The vector field^\ 

dL 



[-ir k S/Hi,L] = [ir n VHi,L] with Hi 



trL 



i+1 



i + 1 



X2) 



(1.3) 



or written out a 
dU 



all commute. 



V)^ - J {IS 



(V) + + J(V) + ) + - ((L*) + J{L%) , L 



(1.4) 



6 Isospectral deformations on tridiagonal symmetric matrices. 

7 V denotes the gradient in the Lie algebra T> with respect to the natural pairing T> 
with V* = V, i.e., VU % = L\ 

8 J{L l ) + means J ((L l ) + ) ; the same for J(L%. 
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The proof of this theorem hinges on a number of Lemmas and the 
Adler-Kostant-Symes theorem 0: 

Lemma 1.2 J is a Lie algebra isomorphism and an order 2 involu- 
tion on T>: 

J 2 = I, (1.5) 
with N+ and M- being the ±1 eigenspaces of J ': 

M± = {a such that Ja = ±a} = {b±Jb with b G V}. (1.6) 

Then T> is decomposed into the Lie algebra A/+ and the (symmetric) 
vector space M- 

V = A/+ + N- , with [N±,N±\ C N+ and [N+,N-\ C A/1, (1.7) 
with 

N± = T>f® (Af±\V ), (1.8) 

and 

A/±\2?o = {a ± suc/i i/iaf a € X>_ or T> + } 

Vf = Vq n A/± = {a ± J a such that a G V }. (1.9) 

Proof : The map J is a Lie algebra homomorphism, such that 
J 2 = I; indeed, upon using J 2 = —I, 

J([a, b}) = J[a, b} T J = -J[a T , b T }J = [Ja 1 J, Jb T J] = [Ja, Jb]. 

In particular, J(a) = ±a and J{b) = ±6 imply 

J[a,b] = [±a,±b] = [a,b], 

leading to the inclusions (1.7); so A/+ is a Lie subalgebra of P, namely 
an infinite-dimensional version of the symplectic algebra sp(oo). The 
second description (1.6) of J\f± follows from 

J{b±Jb) = Jb± J 2 b = ±(b±Jb), for be V; 

and the vector space decomposition V = A/+ + A/1: 

a = i(a + Ja) + ^-(a - Ja). 
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Since JV £ ,V £ J C T> £ ,T>J = T>^ £ , for £ = +,—, 0, we have 

J : V £ — ► P_ £ , 
and so, if a G X>+, we have 6 := J a G £>_ and 

A/± D a± J a = J 2 a ± J a = ±(b ± Jb) with b=Ja£V_. 
Therefore 

{6 ± Jb such that b G T>-} = {a ± J a such that 6 G £>+} = 7V ± \X>o 
establishing the first relation (1.9). As to the second, one checks that 

' a - d 2b 

a b \±( 1 \f a b \ f 2c ~(a-d) 



c d I -1 c d \ -1 



(a + d)J 
(1.10) 



Lemma 1.3 The decomposition 

V = ken, 

k = D_ © Dq", n = A/+ = {a swc/i i/iai J7"a = a} = sp(oo), 

is a vector space splitting of T> into Lie subalgebras. Any element 
a G T> decomposes uniquely into its projections onto k and n: 

a = n^a + 7r n a 

= | (a- - Ja+) + i (a - Ja ) j + |(a+ + Ja+) + ^ (a + c7a )| 

Proof : Since [£>-,£>_] C P_, since V-V , T> V- C P_ and 
[V ,V ] = 0, we have that [k,k] = P_ C k, so that k = £>_ + 
and n = A/+ are both Lie subalgebras. 

But also k n n = 0; indeed, a typical element of A/+ has the form, 

(b-+J(b-)) + (d +J(d )) g (M + \v ) + (M+nv ) = (N+\v )+v+, 
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with 6_ £ do € T>q. If this element € k = £>_ + T> , then 

= and thus b_ = 0, but then, since d + J d Q E T>q n T>q = 0. 
Finally, any element a £ T> can be decomposed intc0 

a = a_ + ciQ + a + 



+ |(a+ + ^a+) + ^ (a + Ja ) J 

€ (x>_ + e ((JV+V^o) +^o") = k © n - 



We now state the R-matrix version of the AKS-theorem (see @ 
and |l^]). The substance of the R-matrix extension is that initial 
conditions £(0) can be taken in all of g, instead of merely in k. 

Proposition 1.4 (Adler-Kostant-Symes, Russian flavor) Letg = 
k + n be a (vector space) direct sum of a Lie algebra g in terms of 
Lie subalgebras k and n, with g paired with itself via a non- degenerate 
ad-invariant inner producf**] ( , ); this in turn induces a decomposition 
g = k -1 " + n -1 - and isomorphisms g ~ g*, k -1- ~ n* ; ~ k*. It also 
leads to a second Lie algebra gR ~ g R derived from g, namely: 

gR ■ [x,y]n = ^[Rx,y] + ^[x,Ry] = [vr k x,7r k y] - [vr n x,7r n y], 

with R = vrk — vr n . Ad* ~ Ad-invariant functions ip on g* ~ g Poisson 
commute for the Kostant-Kirillov Poisson structure 

{/, h} R (0 := <£, [Vh(0, VfiOU) on g* R ~ g R . 

The associated Hamiltonian flows are expressed in terms of the Lax 
pai^\ 

1= [-7r k Vp(£),£] = [7r n V^(£),C] for £ £ gR. 
The systems has a solution expressible in two different waysf^\: 

£(i) = Ad xw £ = Adjv-i W £o (1.11) 



'i 



a± = TTv±a, a — 7Tp a. 
g X;Y) = (X,Ad g . 



10 (Ad g X;Y) = (X,Ad g . 1 Y),g£G, and thus ([z,x],y) = (x,-[z,y\). 



11 V(/? is defined as the element in g* such that d<p(£) = (V<p, d£), £ £ g. 
12 naively written Ad K(t) £ = K(t)£ K(t)-\ Ad N -i& = ^- x (t)£ iV(t). 
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witf^ 

K(t) = n Sk e tV ^°\ and N(t) = TTg a e tV ^°\ 

Proof of Theorem 1.1: Identifying k and n of Lemma 1.3, with those of 
proposition 1.4 and unraveling the projection ir^, according to Lemma 
1.3, imply that the vector fields (1.3) and so (1.4) all commute. ■ 



2 The vector fields dm/dtk = A k m + 
mA Tk 

The main claim of this section can be summarized in the following 
statement: 

Theorem 2.1 Consider the skew- symmetric solution 



m - 



) (t) = e E^ moo(0)e E' 



to the commuting equations 
dt k 



A K m 00 + m 00 A Tk , 



(2.1) 



with skew- symmetric initial condition m(0) and its "skew-Borel de- 
composition" 



moo = Q 1 JQ 1T , with Q G Q^. 



(2.2) 



Then the matrix Q evolves according to the equations 

-7T k (QA^ 1 ] 



Oti 



(2.3) 



and the matrix L := QAQ 1 provides a solution to the Lax pair 



8L 

dU 



n n L\L 



(2.4) 



Conversely, if Q €E Gk satisfies (2.3), then moo, defined by (2.2), sat- 
isfies (2.1). 



13 with regard to the group factorization A = irg k A Trg n A. 
14 Remember C/k denotes the group of invertible elements in k. 
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The proof of this theorem hinges on the following proposition: 
Proposition 2.2 For the matrices 

L := QAQ~ l and m := Q _1 JQ~ 1T , with Q € Qk, 
the following three statements are equivalent 

(ii) U + ^CT 1 € n 

(Hi) ^ = A l m + mA T \ 

Whenever the vector fields on Q or m satisfy (i), (ii) or (Hi), then the 
matrix L = QAQ^ 1 is a solution of the Lax pair 

-gj- = [-KkL\L] = TT n L\L 

Proof : Written out, proposition 2.2 amounts to showing the equiv- 
alence of the three formulas: 



(I) |f Q 1 + {(L% ~ J{L\) T J) + \ ((L*) - J({V)o) T j) = 

(II) ( i . + g Q - 1 )_ J ( i . + f<r>) T ^o 



(III) A*m + mA Ti - — = 0. 
V ' dti 

The point is to show that 



(I) + = 0, (I)_ = (II)_ = -J (II + ) T J, (I) = \ (II) , 

Q-\lI)JQ- lT = (III). (2.5) 
The fact that Q £ G k C k = V_ ®V = V_ ® (V nN~) amounts to 

Q + = and (Q + JQ T J) = 0. 

The latter statement implies that Q is lower-triangular and, along the 
diagonal, is composed of 2 x 2 blocks proportional to the identity. 
Hence 

Qo, Qo e5 "> and so (QQ _1 )o g £>5~ 
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and thus J commutes with (QQ 1 )o, yielding: 

J{QQ~ x )oJ = JHQQ' 1 ^ = -(QQ _1 )o- (2.6) 

Also notice that for any matrix A 

{JAJ) ± = JA±J and (JAJ) = JA J. 
At first, one observes that 

(I)+ = 0. (2.7) 

Also 

(II)- = [U + QQ- 1 - J{L l + QQ- l ) T j)_ 

= (D + QQ- 1 - J((L*) + ) T j)_, using (QQ- 1 )+ = 
= (I)- (2-8) 



(II)o = [V + QQ- 1 - jiV + QQ- 1 ) 1 J 



and 



( L 0o- J ( LlT )o J+ (^" 1 )o- J (W" 1 ) T )o J 
(L l ) - J(L lT ) J + 2(QQ~ 1 )o, using (2.6) 

2(I)o (2.9) 



(ii)+ = (l' + qq- 1 - j(v + qq- 1 ) t j) + 

= (L i ) + -J((L%) T J-J((QQ- 1 )-) T J 

= (((L i )+) T -J(L i )_J-J(QQ- 1 )_j) T 

= J (J((L%) T J - J 2 (L%J 2 - J\QQ- 1 )_J 2 ) T J 

= -j((L l + QQ- 1 - J(L\) T J)-) T J 

= -J((I)-) T J. (2.10) 
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Using the definitions L = QAQ' 1 , m = Q~ 1 JQ~ 1T and J 2 = —I, one 
finds 

(III) = AWmA Tl -^ 

= K l m + mA T ' + (Q~ l Q)m + m(Q' 1 Q) J , setting m = Q' 1 JQ~ lT 
= Q-\QVQ- 1 )JQ- lT -Q- 1 (JQ- lT A Ti Q T J)JQ- lT 
+ Q- 1 (QQ- 1 )JQ- lT - Q- 1 (jQ- lT Q T j) JQ~ lT 

= Q^iyV + QQ- 1 -J{V + QQ- 1 ) T j)jQ- lT 

= Q-\ll)JQ- lT . (2.11) 

The five facts (2.7) up to (2.11) imply (2.5) and so the desired equi valences. ■ 



3 Skew-orthogonal polynomials 

As before, consider the skew-symmetric solution 

moo (t) = e^ Ak moo (0)e^ tkATk 

to the commuting equations (2.1). We now perform the "Borel de- 
composition" of : 

moo = Q- 1 JQ T -\ (3.1) 

where Q € Gk (lower-triangular invertible matrices in k). To carry 
out this factorization, we found the following useful algorithm, based 
on the observation that the skew-Borel decomposition is tantamount 
to finding skew-orthogonal polynomials, in the same way that Borel 
decomposing a Hankel matrix is tantamount to finding orthogonal 
polynomials. 

Define the skew-inner product (/, g) between formal Taylor series, 
induced by 

(y\z 3 ) := mj(t) = {moo{t))ij. 
Monic polynomials p^ are said to be skew- orthogonal, when 

({Pi,Pj))o<i,j <oo = Jh, with h € Vq. (3.2) 
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They, in turn, lead to skew-orthonormal polynomials qk(z), 

1 1 

q2n(t,z)= j= — P2n(t,z) and q 2n +i (M) = j= — p 2n +i (t,z), (3.3) 

V ^2n V h,2n 

thus satisfying 

({li,Qj))o<i,j«x> = J- (3.4) 

Indeed, finding a set of skew-orthonormal polynomials is tanta- 
mount to the "skew-Borel decomposition!' (3.1). Using the lower- 
triangular matrix Q of coefficients of the polynomials q{z) = Qxi z ), 
we have 



J = ((m))^ = (((Qxh (QxWijzo = Q (^' z ')) M >o gT = ^ m °^ T > 

(3.5) 

from which (3.1) follows, and vice- versa. 
Given a skew-symmetric matrix 

A = (aij)o<i,j<n-l, even n, 

we define the Pfaffian by means of the formula^ 
pf(A)dx A ... A dx n _i 

n/2 



2 n / 2 (n/2)\ {^2 £ ^ ai o,h a i2,i3--- a 



j dx A ... A dx n -i, 
(3.6) 



so that pf(A) 2 = det A, but with a sign specification for pf(A). Hence- 
forth, we shall be using the notation 

pf(io,ii,...,ie-i), < i < h < ■■■ < ie-i, £ even, 

to denote the pfaffian of the skew-symmetric matrix, formed by the 
io,ii, ...,i£-ith rows and columns of m^. 



15 In the formula below (ig, i\, i n -2, in-i) — c(Q, 1, n— 1), where a is a permutation 
and e(cr) its parity. 
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Theorem 3.1 Consider a semi-infinite skew- symmetric matrix m^, 
satisfying (2.1); then the monic polynomials p^ 

( o 
-Moi 



P2n(t,z) = 



1 



T2n{t) 



Pf 



A»oi 




V -i 



M0,2n 1 \ 
Ml,2n ^ 







v 2n 



v 2n 



P2n+l(t,z) = 



T2n 



Pf 





-Moi 



Mi 




1 M0,2n+1 \ 

2 Ml,2n+1 







y 2n+l 



\— M0,2n+1 — /^l,2n+l 
9 



,2n+l / 



f 2n (t) K Z+ Wl )T 2n (t)p 2n (t, Z ). 



(3.7) 



/orm a skew- orthogonal sequence, (pi,Pj)o<i,j<oo = Jh, where 



h = (ho, ho, h 2 , h 2 , ■ ■■) G T> , with h 2n 



T 2n +2 
T 2n 



f 2n =pf(0,...,2n-l). 



(3.8) 



whereas the qk, defined in (3.3) from the p's, form an orthonormal 
sequence. The matrix Q defined by q(z) = Qx( z ) is the unique solution 
(modulo signs) to the skew-Borel decomposition of m^: 

m^t) = Q~ l JQ T ~ l , with Q € k. 

The matrix L = QAQ^ 1 , also defined by 

zq(t,z) = Lq(t,z), 

and the diagonal matrix h satisfy the equations 

dL 

dti ~~ 

Remark : It is of interest to mention that, except for the second 
expression for p 2n +i in (3.7), the first part of theorem 3.1 is actu- 
ally t-independent. Namely, (3.7) gives a geometric map, but also 
an effective algorithm, to perform the skew-Borel decomposition of a 
skew-symmetric matrix m^, by applying formulae (3.7). 



-ir k L l ,L 



dh 

and h~ l — = 27r k (L 4 ) . 
dti 



(3.9) 



15 



H - v M: Skew-Orthogonal oct 28, 1998 #1 §3, p 



Corollary 3.2 The polynomials (po(t, z),pi(t, z), ...) have the follow- 
ing explicit form in terms of Pfaffians of the moment matrix m and 
the Pfaffian f -functions: 

^ k pf(0,...,k,...,2n) 

P2n{z) = }^ (~ Z > 



0<fc<2n 



p/(0,...,2n-l) 



P2n+l{ ) " \<LJ ] 2 " ' ; " 

(3.10) 

Remark : Taking into account the identity pf '(odd set) = and the 
convention 

p/(0,...,fc,...,2n,2n + l) = -p/(0,...,2n-l) for fc = 2n+l, 
the polynomials can be written as follows: 

p 2n {t,z) = — - (-z) k pf(0,...,k,...,2n,2rf+l) 
P2n+i(t,z) = ^— Yl (-z) k pf(0,...,k,...,2n,2n + l). 

T2n[ l> 0<fe<2n+l 

(3.H) 

The proof of this theorem hinges on a number of lemmas, involving 
properties of Pfaffians: 

Lemma 3.3 Consider an arbitrary skew- symmetric matrix a = (aij)i,j>c 
For odd t and fixed i > i, the following holds: 

]T (-l) k a ki pf(0,...,k,...,e-l) =pf(0,...,£-l,i). 

o<k<e-i 

Proof : On the one hand 

det(a^) <^ = pf 2 (0, 1, ...,£) (3.12) 
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and, on the other hand, using 



det 



/ 



— «0r 







= p/(0,...,r, k)pf(0,...,r,£), 



* / 



and expanding according to the last row and column 

det(aij)o<ij<e 

E (-iy +i ae t a ] epf(0,..J,...,e-l) P f(0,...,j,...,e-l) 

0<i,j<£-l 

]T (-l)W/(0,...,i,...,€-l) 

,0<i<£-l 

( E (-1)^/(0, 
\0<?<*-i 

]T (-l)W/(0,...,i,...,^-l)] . (3.13) 

V 0<i<£-1 / 

Comparing the two expressions (3.12) and (3.13) and taking into ac- 
count the sign, lead to the claim of lemma 3.3. ■ 

Remark : Lemma 3.3 will be applied as follows (for odd £): 

E (- 1 ) ka ik,ilPf( i 0, ■■; h, ■■; U-l) = f>/(«0, -, U-l,h)- 

o<k<e~i 



Lemma 3.4 The polynomials Pk{z) defined by the determinantal ex- 
pression (3.7) equal the expressions (3.10) of Corollary 3.2 and satisfy 

<JPi(z), Pj {z)) itj > = Jh (3.14) 

Proof : At first notice that the leading coefficients of the pfaffians 
in the polynomials (3.7) equal pf(0, 1, 2n — 1) = T2 n (t). Then apply 
Lemma 3.3 (remark following the lemma) to the first case, setting the 
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a's equal to /it's from to 2n and a^n+i = (1 — &k,2n+i)z k - Also apply 
Lemma 3.3 to the second case, by setting the a's equal to fi, from 
to 2n + 1, skipping 2n, and a^^n = (1 — <>k,2n)z k . Then, we find: 

f2n(t)P2n(t, z) = ]T (-z) k pf(0, .... k, 2n - 1, 2n, 2?Tf 1) 

0<fc<2n 

= X! (-^)V(0,...,fc,...,2n) 

0<fc<2n 

T2n(t)p2n+i(t,z) = Yl (-z) k pf(0,...,k,...,2n-l,2n,2n + l) 

0<fc<2n+l 

E (-«) fc p/(0,...,fc,...,2^,2n + l) 

0<fc<2n-l 

+ z 2n+ip/(0,...,2n- 1). 

(3.15) 

Using the expressions (3.15) and Lemma 3.3, one computes the inner- 
product 

/ r \ i\ < i\fc p/(0, fe, 2n) 

(n,W " ) = „£„ ( " 1)w V/(o..^»-D 

p/(0,...,2n,Q 
p/(0,...,2n-l)' 

which vanishes for < z < 2n, and so 

(P2n(^),^) =0 for < i < 2n; 

therefore 

(P2n(*),Pi(*))=0 for 0<i<2n. 
Moreover, using the expression for P2n+i and Lemma 3.3, 

j ^ fc p/(0,...,fc,...,2n,...,2n + l) 
< ^ lW '' > = ^^o^J P/(0,..,2n-l) 

p/(0,...,2^,2n+M) 
p/(0,..,2n-l) 
for < i < 2n+ l,i ^ 2n 



= < 



p/(0,...,2ra,2ra+l,2n) _ p(0,...,2w,2ra+l) _ f 2 n+2(t) _ _Z 
P /(0,...,2n-l) ~ p/(0,...,2n-l) ~ f 2 „(t) ~ 

for i = 2n. 
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and so 

(p2n+i(z),Pi(z)) = S it 2nh2n for < i < 2n + 1; 
the skew-symmetry of (, ) does the rest. ■ 

Lemma 3.5 We have 

d 1 

— pf(i ,...,ie) = ^2pf(i ,...,i r + k,...,i e ). 

k r =0 



Proof : The differential equations (2.1) for m^t) read in terms of 
the elements pif. 

Fill : 

+ Vij+k', (3.16) 



dpij 



dt k 

hence, setting a(i) = i', we compute 

— pf(i ,...,i n -i) 
dtk 

t/i/Zj.Z a re {0,...,„-2} 



n-l 

= ^pf{io,-.-,ir + k,...,i n -i), 

r=0 

upon reorganizing the first sum on the right hand side appropriately^ 
Lemma 3.6 The polynomials p k of Theorem 3.1 satisfy: 

f^(Jr +z ) f2 " (t)P2 " (z)=P2 " +i(z) - 

Proof : Using the explicit form (3.10) of p2 n and Lemma 3.5, one 
computes 

d 



2n 

-E(-^) fc+ W(0,l,...,^,...,2n) 



k=0 
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2n 

• £(-^) fc (Rf (°> k - 1, 2n) + p/(0, fc, 2n, 2n + 1)) 

fc=0 



= z 2n+1 p/(0, 1, 2n - 1) + z 2n (-p/(0, 2n^ 1, 2n) +p/(0, 2r^ 1, 2n)) 

2n-l 

+ ^ (-z) k pf(0, k, 2n, 2n + 1) + .., 2n, 2n + 1) 
l 

= T 2 „(t)p2n+l(^), 

ending the proof. ■ 

Lemma 3.7 The decomposition 

moo = Q- 1 JQ T -\ with Q e Q k , (3.17) 

is unique (modulo signs) and is completely given by the skew-orthonormal 
polynomials qk(z). 

Proof : In view of relations (3.5), the skew-orthonormal polyno- 
mials (3.3) provide a solution to the Borel decomposition = 
Q~ 1 JQ T ~ 1 , with Q £ g^. We now show this Q 6 5k is unique! 
To do so, we proceed by induction on g% = (Qx( z ))k- Assuming the 
existence of another matrix Q E Q^, such that rrioo = Q~ 1 JQ T ~ 1 , 
we must show q = q. Indeed moo = Q~ 1 JQ T ~ 1 leads unambiguously 
(except for the sign) to 

Qo = Hm 12 , and q x = fi i /2 z. 

But, for some a% and (3i, we would have, taking into account the 
absence of the z 2n -term in fa n +i, 

2n 2n-l 

fan = ^oiiqi and q 2n +i = a 2n q2n+i + ^ ( 3 - 18 ) 


Then using the inductive step, namely fa = q^ for < k < 2n — 1, 
and the skew-orthonormality, we compute for £ < n, 

= {fan, fat) = (<?2n,<?2<?) = -«2<?+l 

= {fan,fat+l) = {fan,<lU+l) = <*2l 

= {fan+l,fal) = {fan+l,<l2l) = ~@2t+l 

= (?2ti+1) 52^+1 ) = (<?2n+l, 92^+1 ) = thl, 
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and so, since also /?2n = 0, 

hn = «2n<?2n and q 2 n+l = ^n^n+l- 

Finally 

1 = (hmhn+i) = al n {q 2n ,q2n+i) = a| n , 

leading to a\ n = 1. ■ 

Proof of Theorem 3.1 and Corollay 3.2 : The proof follows at once 
from lemmas 3.3 up to 3.7. The fact that Q and L satisfy the differ- 
ential equations (2.3) and (2.4) follows from the fact that the matrix 
moo satisfies the differential equations 

dm/dt k = A k m + mA Tk , 

according to Theorem 2.1. 

The second relation in (3.9) is established as follows: defining the 
semi- infinite lower-triangular matrix P (with l's along the diagonal) 
by means of the monic skew-orthogonal polynomials p(z) = Px{ z ) 
and differentiating P = h l / 2 Q with regard to t^, one obtains 

h l ' 2 QQ- l h- 1 ' 2 = -\{\ogh)- + PP~\ 

Then, using proposition 2.2, (formula (II)), one finds: 

= ^/ 2 (II) ^ 1/2 = h 1/2 ((L k ) - J(L kT ) J + 2(QQ' 1 ) )h- 1 / 2 
= h 1 ' 2 !^)^- 1 ' 2 - (log ft)' + 2(P P-^o 
= 2^ k (L fc ) - (log/*)' , 

using the commutation of 7Tk(£ fc )o with h 1 / 2 G T>q , together with 
(P P- 1 ^ = 0. ■ 
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